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(57) Abstract 

An adaptive reception 
and transmission technique 
according to one embodiment 
of the invention offers the 
advantages of adaptive 
reception transmission 
using feedback without 
the associated mobile 
radio complexity increase 
and information capacity 
penalty. The technique has 
been developed to exploit 
structured variation which 
occurs in the multipath 
fading present in the wireless 
antenna array channel. Thus, 
multipath propagation effects 
are explicitly accounted for 
in the problem approach. 
The technique is blind in that 
the antenna beam is formed 
in the absence of explicit 
knowledge of the array 
geometry, and without the 
necessity of array calibration 
or mobile feedback. The 
basic approach is to estimate 
to optimum receive and 

transmit antenna beam pattern based on certain statistical properties of the received antenna array signals. The optimum receive and 
transmit beam pattern is found by solving an optimization equation. The adaptive transceiver system is suitable for use in conjunction 
with either a diplexed transmit/receive antenna array, or with separate transmit and receive arrays. 
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WIRELESS COMMUNICATION NETWORK USING TIME-VARYING 
VECTOR CHANNEL EQUALIZATION FOR ADAPTIVE SPATIAL 

EQUALIZATION 

RELATED APPLICATIONS 
5 This application claims priority to the following applications and incorporates 

these applications by reference: 

U.S. Ser. No. 08/394,652 filed on February 22, 1995; 
U.S. Ser, No. 08/491,044, filed on June 16, 1995; and 
U.S. Prov. App. No. 60/005,647 filed on October 19, 1995. 

10 FIELD 

The present invention relates to a wireless communication network that uses 
time- varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
quality of communication. 

15 BACKGROUND 

Within wireless mobile communication systems, procedures for compensating 
for multipath effects are desirable. Unlike "Iine-of-sight" radio links, a number of 
signal transmission paths typically comprise each wireless communication channel, 
hence the term "multipath." An increase in primary path communication energy is 

20 desired with an attendant reduction in the interference energy radiated to mobile 
users over other non-primary paths. Often this increase in primary path 
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communication energy is achieved through generation of spatially selective, directive 
transmission beam patterns. 

Some systems employ directive antennas at base station sites to increase the 
signal level received from and transmitted to each mobile user relative. Other 

5 systems use an antenna array that is used to form beams to increase the signal level 

received from and transmitted to each mobile user. In the systems with an antenna 
array, often called a phased array, weights are used to account for various angles 
that the beam will provide. Problems occur when severe multipath signals are 
created by a plurality of obstructions such as buildings and mountains. When severe 

10 multipath occurs, many systems suffer severe quality degradation. Known systems 

cannot sufficiently account for severe multipath and the result is reduced quality 
communication. 

Figure 1 shows an illustrative representation of a wireless "multipath" 
communication channel between a base station 2 and a remote mobile user 4. The 

15 various signal transmission paths within each such multipath communication channel 

arise from reflection of the transmitted signal by dominant reflection surfaces 6, and 
by minor reflection surfaces 12, between the base station 2 and remote mobile user 
4. Accordingly, techniques for improving signal reception in line-of-sight radio 
systems are often not directly applicable to multipath signal environments. For 

20 example, in line-of-sight system the "gain" of an antenna typically varies inversely 

to antenna beam width. However, if a given antenna beam width is made less than 
the angular region encompassing the various signal paths comprising a multipath 
communication channel, further reduction in the beam width may reduce the energy 
radiated along some of the angular paths. In turn, this may actually decrease the 

25 effective time average gain of the antenna. 

These systems that perform blind transmission fail to employ a sufficient 
method for determining the received angle of the mobile station signal. As a result, 
known transmission beamforming techniques often use a trial and error approach 
that works barely satisfactorily, and does not always provide high quality 

30 communication. 
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RELATED ART 

Within wireless mobile communication systems, three techniques have been 
developed for improving communication link performance using directive transmit 
antennas: (i) selection of a particular fixed beam from an available set of fixed 
5 beams, (ii) adaptive beam forming based on receive signal angle estimates, (iii) 

adaptive transmission based on feedback provided by the remote mobile user, and 
(iv) adaptive transmit beam forming based upon the instantaneous receive beam 
pattern. Each of these techniques is described briefly below. 

In the first technique, one of several fixed base station antenna beam patterns 

10 is selected to provide a fixed beam steered in a particular direction. The fixed 
antenna beams are often of equal beam width, and are often uniformly offset in 
boresight angle so as to encompass all desired transmission angles. The antenna 
beam selected for transmission typically corresponds to the beam pattern through 
which the largest signal is received. The fixed beam approach offers the advantage 

15 of simple implementation, but provides no mechanism for reducing the signal 

interference power radiated to remote mobile users within the transmission beam of 
the base station. This arises because of the inability of the traditional fixed beam 
approach to sense the interference power delivered to undesired users. 

The second approach involves "adapting" the beam pattern produced by a 

20 base station phase array in response to changing multipath conditions. In such 
beamforming antenna arrays, or "beamformers", the antenna beam pattern is 
generated so as to maximize signal energy transmitted to ("transmit beamforming"), 
and received from ("receive beamforming"), an intended recipient mobile user. 
While the process of transmit beamforming to a fixed location over a 

25 line-of-sight radio channel may be performed with relative ease, the task of 

transmitting to a mobile user over a time-varying multipath communication channel 
is typically considerably more difficult. One adaptive transmit beamforming 
approach contemplates determining each angle of departure (AOD) at which energy 
is to be transmitted from the base station antenna array to a given remote mobile 

30 user. Each AOD corresponds to one of the signal paths of the multipath channel, 
and is determined by estimating each angle of arrival (AOA) at the base station of 
signal energy from the given user. A transmit beam pattern is then adaptively 
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formed so as to maximize the radiation projected along each desired AOD (i.e, the 
AOD spectrum), while minimizing the radiation projected at all other angles. 
Several well known algorithms (e.g., MUSIC, ESPRIT, and WSF) may be used to 
estimate an AOA spectrum corresponding to a desired AOD spectrum. 

Unfortunately, obtaining accurate estimates of the AOA spectrum for 
communications channels comprised of numerous multipath constituents has proven 
problematic. Resolving the AOA spectrum for multiple co-channel mobile units is 
further complicated if the average signal energy received at the base station from any 
of the mobile units is significantly less than the energy received from other mobile 
units. This is due to the fact that the components of the base station array response 
vector contributed by the lower-energy incident signals are comparatively small, thus 
making it difficult to ascertain the AOA spectrum corresponding to those mobile 
units. Moreover, near field obstructions proximate base station antenna arrays tend 
to corrupt the array calibration process, thereby decreasing the accuracy of the 
estimated AOA spectrum. 

In the third technique mentioned above, feedback information is received at 
the base station from both the desired mobile user, and from mobile users to which 
it is desired to minimize transmission power. This feedback permits the base station 
to "learn" the "optimum" transmit beam pattern, i.e., the beam pattern which 
maximizes transmission to the desired mobile user and minimizes transmission to all 
other users. One disadvantage of the feedback approach is that the mobile radio 
needs to be significantly more complex than would otherwise be required. 
Moreover, the information carrying capacity of each radio channel is reduced as a 
consequence of the bandwidth allocated for transmission of antenna training signals 
and mobile user feedback information. The resultant capacity reduction may be 
significant when the remote mobile users move at a high average velocity, as is the 
case in most cellular telephone systems. 

The fourth conventional technique for improving communication link 
performance involves use of an optimum receive beam pattern as the preferred 
transmission beam pattern. After calibrating for differences between the antenna 
array and electronics used in the transmitter and receiver, it is assumed that the 
instantaneous estimate of the nature of the receive channel is equivalent to that of the 
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+' transmit channel. Unfortunately, multipath propagation and other transient channel 

phenomenon can substantially eliminate any significant equivalence between 
frequency-duplexed transmit and receive channels, or between time-division 
duplexed transmit and receive channels separated by a significant time interval. As 
5 a consequence, communication link performance fails to be improved. 

SUMMARY 

The present invention relates to a wireless communication network that uses 
time- varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
10 quality of communication. 

Accordingly, it is an object of the invention to provide an adaptive receive 
and transmit beamforming technique that improves communication quality by 
utilizing the received signal energy from remote users. 

It is another object of the invention to provide an adaptive receive and 
15 transmit beamforming technique which accounts for the presence of multipath 

fading inherent in the communication channel. 

It is yet another object of the invention that the beamforming technique be 
independent of antenna array geometry, array calibration, or of explicit feedback 
control signals from remote users. 
20 It is another object of the invention to provide adaptive receive and transmit 

beam forming that improves signal quality received by a desired user and while 
simultaneously reducing interference energy received by other undesired users so as 
to, within a cellular communication network, improve communication traffic 
capacity, and/or to increase base station coverage area, and/or to improve call 
25 quality. 

An adaptive reception and transmission technique according to one 
embodiment of the invention offers the advantages of adaptive reception 
transmission using feedback without the associated mobile radio complexity increase 
and information capacity penalty. The technique has been developed to exploit 
30 structured variation which occurs in the multipath fading present in the wireless 

antenna array channel. Thus, multipath propagation effects are explicitly accounted 
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for in the problem approach. The technique is blind in that the antenna beam is 
formed in the absence of explicit knowledge of the array geometry, and without the 
necessity of array calibration or mobile feedback. The basic approach is to estimate 
the optimum receive and transmit antenna beam pattern based on certain statistical 

5 properties of the received antenna array signals. The optimum receive and transmit 
beam pattern is found by solving an optimization equation. The adaptive transceiver 
system is suitable for use in conjunction with either a diplexed transmit/receive 
antenna array, or with separate transmit and receive arrays. 

Advantages of the invention include reduced interference and improved 

10 communication quality. 

BRIEF DESCRIPTION OF THE FIGURES 
Additional advantages of the invention will become apparent upon reading 
the following detailed description and upon reference to the drawings, in which: 

Figure 1 is an illustrative representation of a multipath propagation channel 
15 between a base station and a remote user station; 

Figure 2A is a block diagram of the physical implementation of a 
beamforming network configured to perform adaptive beam forming in accordance 
with the present invention; 

Figure 2B is a block diagram of the physical implementation of a 
20 beamforming network in which a diplexor is employed to allow the antenna array to 
be used for both transmit and receive operation; 

Figure 3 is a functional block diagrammatic representation of a beamforming 
network of the present invention; 

Figure 4 is a graph depicting estimated signal SINR versus adaptation time 
25 for signal 1 ; 

Figure 5 is a graph depicting estimated signal SINR versus adaptation time 
for signal 2; 

Figure 6 is a graph depicting estimated signal SINR versus adaptation time 
for signal 3; and 

30 Figure 7 is a graph showing theoretical performance bounds for variable A. 



-6- 



WO 98/18271 



PCTYUS97/18196 



DETAILED DESCRIPTION 
The present invention relates to a wireless communication network that uses 
time-varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
5 quality of communication. 

The exemplary embodiments are described herein with reference to specific 
configurations and protocols. Those skilled in the art will appreciate that various 
changes and modifications can be made to the exemplary embodiments while 
remaining within the scope of the present invention. 

10 The problem we consider is blind adaptive spatial equalization with a 

multiple antenna base station array for the purpose of increasing reverse link spectral 
efficiency (Erlangs/Hz/base station) or range by canceling interference and 
enhancing the desired signal over noise. In blind equalization, the equalizer weight 
vector adaptation relies on the structure of the transmitted signal rather than a known 

15 training signal. While these techniques may be applied at the portable unit, this 

present embodiments focus on spatial equalization at the base station. 

The communication channel formed between a portable wireless user and a 
wireless base station varies with time. When the portable moves with high velocity, 
the time variation in the channel can become severe. Over the past decade, a large 

20 body of work has been devoted to blind spatial equalization in time invariant 

channels. However, less work has been reported for large Doppler-spread 
conditions. For certain communication signals, including narrow-band FM, the 
known existing equalizer algorithms fail to converge in rapidly fading environments. 
This occurs because the channel changes significantly over the time interval required 

25 to determine the next equalizer weight vector update. 

In most conventional adaptive algorithms, the signal quality criteria is 
minimum mean squared error (MMSE). We will show that this criteria may not be 
the best choice in rapidly fading channels. We will also show that the signal to 
interference plus noise ratio (SINR) criteria can be more robust in rapidly fading 

30 channels. With this information and motivation, an algorithm is developed which 
adaptively optimizes an interesting variation of the conventional SINR criteria. 
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The exemplary equalizer decomposes the received vector signal into a desired 
signal subspace and an interference plus noise subspace. In order to accurately 
estimate these two subspaces, we exploit the inherent time correlation structure that 
exists in the fading channel to obtain instantaneous blind estimates of the desired 
5 channel. Supplement 1 is provided for additional information and includes 

additional background and introduction material. 

I. Overview of Beamforming Network 

Turning now to Figure 2 A, a block diagram is shown of the physical 
organization of a beamforming network 50 configured to perform adaptive beam 

10 forming in accordance with the present invention. In an exemplary embodiment the 

beamforming network 50 is disposed within a base station of a cellular 
communications network, in which is included a transceiver comprised of a radio 
frequency (RF) transmitter 52 and an RF receiver 54. 

In the embodiment of Figure 2A, a base station antenna array 56 serves to 

15 produce independent transmit and receive antenna beams for facilitating 

communication with one or more mobile units (not shown). The term "receive 
channel vector" is employed to indicate that each antenna element within the base 
station antenna array 56 will form a propagation channel to a given remote user. 
The composite array channel may be represented as a vector having elements 

20 corresponding to each individual antenna channel. As is described herein, statistical 
characterization of the receive channel vector provides information which may be 
used by the base station to determine an "optimal" transmit beam pattern, i.e., a 
transmit beam pattern which maximizes the average signal-to-interference power 
delivered to a given mobile user. This obviates the need for the mobile unit to 

25 provide feedback information to the base station relating to propagation 

characteristics of the transmit channel. This in turn simplifies implementation of the 
mobile unit, and preserves channel resources for data transmission by eliminating the 
need for "mobile unit feedback" relating channel characteristics to the base station. 
As is indicated by Figure 2 A, a diplexer 58 can be employed to allow the 

30 antenna array 56 to be used for both transmit and receive operation by isolating the 

RF receiver 54 from the RF transmitter 52. A receive channel beamformer 60 
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cooperates with the RF receiver 54 to adaptively optimize the receive antenna beam 
in order to improve received signal quality. Similarly, a transmit channel 
beamformer 64 cooperates with the RF transmitter 52 to adapt the transmit antenna 
beam to optimize some characteristic of transmission quality. In an exemplary 
5 embodiment the transit channel beamformer 64 and receive channel beamformer 60 
are each implemented as a special purpose digital signal processor (DSP). 

In another embodiment of the invention, distinct antenna arrays are used for 
signal reception and transmission as illustrated in Figure 2B. In the embodiment of 
Figure 2B, a diplexer is not required since a dedicated receive antenna array (not 

10 shown) is connected to the RF receiver 54, and a dedicated transmit antenna array 

(not shown) is connected to the RF transmitter 52. The receive and transmit antenna 
arrays are designed to provide identical radiation characteristics when operated at the 
receive and transmit frequencies, respectively. Accordingly, in many instances the 
physical geometries of the transmit and receive antenna arrays are simply physically 

15 scaled to account for the fractional difference in the receive and transmit RF 

wavelengths. The embodiment of Figure 2B substantially eliminates the potential 
introduction of error arising from use of a single antenna array and diplexer. 

Turning now to Figure 3, a functional block diagrammatic representation is 
provided of a beamforming network of the present invention. In Figure 3, solid 

20 lines are used to represent functional elements and dashed lines are employed to 
identify the physical components of Figure 2. The RF receiver 54 is functionally 
comprised of a low-noise amplifier (LNA) network 110, and a complex 
downconverter 112. The complex downconverter 112 is disposed to frequency 
downconvert the received RF signal energy after processing by the LNA network 

25 110. The downconverted signal energy is then digitally sampled and provided to a 

receive channel weighting module 1 16 of the receive channel beamformer 60. The 
weights applied by the receive channel beamformer 60 to each of the M 
downconverted antenna element outputs x Rm (k), m =1 to M, of the complex 
frequency downconverter 1 12 are determined by a receiver weight vector adaptation 

30 module 118. In the exemplary embodiment the receiver weight vector adaptation 
module 118 determines a receive channel weight vector, w R , which maximizes the 
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signal quality received over the desired inbound frequency channel. The weight 
vector calculation is performed according to that described in Supplement 1. 

In the embodiment of Figure 3, a vector channel covariance estimator 140 
within the transmit beamformer 64 operates to produce a statistical characterization 
of a receive channel vector using: (i) the outputs x Rm (k), m =1 to M, of the 
complex frequency downconverter 112, and (ii) an estimate of the desired signal 
S r (k) generated in the receive channel beamformer 60. For present purposes the 
receive channel vector may be viewed as being representative of the multipath 
communications channel from a mobile user (not shown in Figure 3) to the antenna 
array 56. In an exemplary embodiment the statistical characterization carried out 
within the covariance estimator 140 yields an estimated receive channel covariance 
matrix used during the transmit beamforming process. Throughout the following 
description, scalar quantities are represented using lower case characters of standard 
type, vector quantities are represented using lower case characters of bold type, and 
matrix quantities are represented using upper case characters of bold type. 

Within the transmit channel beamformer 64, an optimal transmit beam 
pattern weight vector is generated by a transmit channel weight vector adaptation 
module 150 based on the results of the statistical characterization of the receive 
channel vector. This weight vector, w T (t), optimizes a particular aspect (e.g., 
average desired signal to undesired interference ratio) of the signal energy within the 
transmission range of the base station array 56. In the exemplary embodiment, the 
optimal transmit beam pattern weight vector is generated using the estimated desired 
receive channel covariance matrix, R^k), and undesired interference covariance 
matrix, R^k), both of which are compiled within the covariance estimator 140. 

As is indicated by Figure 3, the signal information (S T ) to be transmitted to 
the desired mobile radio unit is used to modulate a digital baseband carrier within a 
modulator 154. The modulated signal is then applied to a transmit channel 
weighting module 158 disposed to weight, on the basis of the optimized transmit 
pattern weight vector, the input signals corresponding to each element of the 
antenna array 56. The weighted set of input signals produced by the weighting 
module 158 are then upconverted in frequency within a complex frequency 
upconverter 160 of the RF transmitter 52. The resultant frequency-upcon verted 
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signals are then amplified by a bank of power amplifiers 164, and provided to the 
antenna array 56 for transmission via diplexer 58. 

In the exemplary embodiment an improved estimate of the received signal is 
obtained through utilization of a demodulator 170 and remodulator 172. The 
received signal is demodulated within demodulator 170 in order to recover the 
essential characteristics of the modulating signal. In the case of an analog FM 
signal, this involves recovery of the FM waveform. In the case of a digitally 
modulated signal (BPSK, FSK, QPSK, etc.), the demodulator 170 forms hard 
decisions as to the value of each digital symbol. The demodulated signal is then 
processed based upon some predefined characteristic of the signal and modulation 
pattern. For example, a demodulated analog FM signal could be lowpass filtered 
based upon a known signal information bandwidth as a means of obtaining an 
improved post-demodulation signal estimate. In the case of digital modulation, error 
correction could be implemented in order to remove bit errors, thus improving 
estimated signal quality. In addition, training signals (i.e., pilot tones, SAT tones, 
etc.) may optionally be employed in lieu of, or in conjunction with, the 
aforementioned "blind" techniques. 

Again referring to Figure 3, the processed demodulated signal is then used by 
demodulator 172 to remodulate an RF carrier, thereby producing an improved 
modulated signal estimate. The improved signal estimate is then used by the 
receiver weight vector adaptation block 118 and the covariance estimate 140. Other 
techniques, which do not rely upon such a demodulation/remodulation procedure, 
can be devised for obtaining a sufficiently accurate received signal estimate. Figure 
3 simply illustrates a particular exemplary embodiment incorporating a demodulator 
and remodulator. 

In the present embodiment, the demodulator 170 and remodulator 172 or the 
receive channel beamformer 60 are operative to produce a received signal estimate 
s r . The quantity S r is then employed by the covariance estimator 140 to estimate the 
covariance matrix of the receive channel. The receive channel beamformer 60 could 
of course be replaced by numerous alternative structures including, for example, 
multi-antenna sequence estimators, maximum likelihood receiver structures, multiple 
element decision feedback equalizers, or similar alternative described in Supplement 
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1 . Any of these alternative structures could also be used to provide the quantity s r 
for use in estimating the received channel covariance statistics. 

II. Blind Beamforming Using Multipart Signals 

As mentioned above, an initial step within the present invention involves 
5 characterizing a receive channel vector representative of the multipath 

communications channel from a mobile user to a base station antenna array 56. In 
one embodiment, this is accomplished as described in Supplement 1 . This 
embodiment employs a five step approach: 

1) estimate the present signal sequence by passing the equalizer output 
10 through a decision device; 

2) estimate the present channel state with an optimal estimation filter; 

3) estimate the received interference plus noise signal sequence and use this 
error signal to update the interference covariance estimate; 

4) update the spatial equalizer weight vector given the estimated desired 
15 channel and interference covariance; and 

5) repeat the process with new desired signal estimates. 

III. Alternate Embodiments 

Although the embodiments of the invention presented herein have been 
described in terms of optimization from a base station perspective, the teachings of 

20 the invention are equally applicable to optimum formation of the mobile unit antenna 
pattern. That is, the mobile unit antenna pattern is formed such that a desired level 
of power is transmitted to the one or more base stations with which the mobile unit 
is in communication, and so that transmission power is minimized to the other base 
stations within the communication network. Again, an initial step in this antenna 

25 pattern optimization procedure involves statistically characterizing the mobile unit 
receive channel vector. The results of this statistical characterization are then used 
to determine the beam pattern weight vector which maximizes a predetermined 
quality parameter (e.g., signal to noise ratio) of the signal received by the intended 
base station(s), while minimizing signal transmission to other base stations. At the 
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mobile station, the transmit antenna beam pattern is then generated in accordance 
with the transmit beam pattern weight vector. 

Accordingly, in the following claims the term "central communications 
station" is intended to refer to either the base station or mobile unit configured to 
5 generate an optimized antenna beam pattern in accordance with the invention. When 
a base station comprises the central communications station, the mobile unit(s) in 
communication therewith are referred to as "remote communications station(s)". 
Conversely, when a mobile unit comprises the central communications station, the 
receiving base station(s) are identified as remote communications station(s). 
10 While the present invention has been described with reference to a few 

specific embodiments, the description is illustrative of the invention and is not to be 
construed as limiting the invention. Various modifications may occur to those 
skilled in the art while remaining within the true spirit and scope of the invention as 
defined by the appended claims. 

15 SUPPLEMENT 1 

New blind adaptive spatial equalization techniques are proposed for severe 
multipath fading with co-channel interference. In rapidly fading environments, 
equalizer performance suffers with existing algorithms because the desired and 
interference channels change significantly over the update time interval. It is shown 

20 that for fast fading channels, conventional squared error equalizer cost functions can 

yield poor results. A new algorithm is proposed based on a variation of the SINR 
cost function coupled with optimal time varying channel estimation to improve 
signal tracking. In simulations, the algorithm achieves theoretically optimum time 
varying equalizer performance in extreme interference channels with narrow band 

25 FM signals and high velocity portables. These results are contrary to the widespread 
belief that high performance blind adaptive array processing is not possible for 
AMPS with severe combinations of Doppler spread, angle spread, and interference 
power. 
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1. INTRODUCTION 
The problem we consider is blind adaptive spatial equalization with a 
multiple antenna base station array for the purpose of increasing reverse link spectral 
efficiency (Erlangs/Hz/base station) by canceling interference and enhancing the 
5 desired signal over noise. In blind equalization, the equalizer weight vector 

adaptation relies on the structure of the transmitted signal rather than a known 
training signal. While these techniques may be applied at the portable unit, this 
paper focuses on spatial equalization at the base station receiver. 

The communication channel formed between a portable wireless user and a 

10 base station varies with time. When the portable moves with high velocity, the time 

variation in the channel can become severe. Over the past decade, a large body of 
work has been devoted to blind spatial equalization in time invariant channels. 
However, less work has been reported for large Doppler-spread conditions. For 
certain communication signals, including narrow-band FM, the known existing 

15 equalizer algorithms fail to converge in rapidly fading environments. This occurs 

because the channel changes significantly over the time interval required to 
determine the next equalizer weight vector update. 

In most conventional adaptive algorithms, the signal quality criteria is 
minimum mean squared error (MMSE). We will show that this criteria may not be 

20 the best choice in rapidly fading channels. We will also show that the signal to 
interference plus noise ratio (SINR) criteria can be more robust in rapidly fading 
channels. With this motivation, an algorithm is developed which adaptively 
optimizes an interesting variation of the conventional SINR criteria. The proposed 
equalizer decomposes the received vector signal into a desired signal subspace and 

25 an interference plus noise subspace. In order to accurately estimate these two 

subspaces, we exploit the inherent time correlation structure that exists in the fading 
channel to obtain instantaneous blind estimates of the desired channel. 

The multiple antenna wireless channel model is briefly reviewed in Section 2. 
The spatial equalization problem is developed and a solution is proposed in Section 

30 3. Simulation results demonstrating the performance of the adaptive spatial 
equalizer with multiple co-channel narrow-band analog FM signals in a high 
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Doppler spread, high angle spread fading environment are presented in Section 4. 
Conclusions are offered in Section 5. 



The wireless communication channel is characterized by scattering and 
reflecting bodies which cause severe multipath. This is illustrated in Figure 1 . The 
channel is also time-varying when the receiver is in motion. The time-variation is 
characterized by rapid amplitude fluctuations and therefore is called a fading 
channel. These short-term complex amplitude variations are often well modeled by 
a multiplicative Gaussian random process. The channel also exhibits long-term 
variations which are due to shadowing by large physical features. This variation is 
much less structured and is sometimes well modeled as a log-normal multiplicative 
random variable. 

Let s 0 (t) be the complex analytic baseband representation of the transmitted 
signal. If the remote transmitter has a single antenna and the wireless base station 
receiver has multiple antennas, then the direction of arrival of the 1 th path form the 
receiver (given by 6 01 ) is of importance. The baseband vector output of the base 
station receiver dues to the 1 th propagation path is 



where a(6 0 ^ denotes the base station array response to the path arriving at the angle 
6 0I . The composite complex-valued propagation path strength p 0 ,i(t) is given by 



where a 0 ,(t) is the complex fast fading parameter, T, is the shadow fading loss 
factor, d 0 is the distance from the mobile to the base station and r\ is the path power- 
loss exponent. In Equation (1) an inherent assumption is that each propagation path 
amplitude is frequency invariant over the bandwidth of the signal and that the 
relative time delay between different arriving paths is small compared to the inverse 
signal bandwidth. Below, it is further assumed that the array response a(8 0 J) is 



2. CHANNEL MODELS 



x J (t)=a(0 OtJ )P o-J (t)s o (t) 




1 



(2) 
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invariant over the signal bandwidth. This set of conditions is commonly referred to 
as the narrow band assumption. 

The received time-domain signal vector is the sum of all dominant reflector 
paths and the additive receiver noise n(t) 

xit) ^i^Xjit) + n(t) =£a<0 o#J )P Oti (t)s o (t) +/)(t) < 3) 



Clearly from the above model, the time- varying response of the continuous-time 
desired portable channel is 

v t > 5 4 a <© 0 .i>Po.j (t) (4> 



The channel given in Equation (3) is sampled at a rate of T seconds. The 
output samples can be rewritten as, 

x * =h o,* s o,* + /3 * (5) 

with the discrete time channel response given by 

10 In the presence of multiple user co-channel interference, we have the received signal 
vector 

When portables move at high velocities, giving rise to large Doppler 
frequency shifts, the channel vector complex M-space direction and magnitude 
become rapidly varying functions of time. The propagation environment determines 
15 the exact nature of the vector variation. The central path angle and relative spread 

in path angles (angle spread) for each channel determines the Eigenvectors and 
Eigenvalues of the time averaged vector subspace spanned by the channel. The 
Doppler spread determines the rate of time variation within the channel subspace. 
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3. SPATIAL CHANNEL EQUALIZATION 
The techniques proposed in this paper do not consider conventional antenna 
concepts such as beam width, gain and side lobe rejection. Instead, the antenna 
array is viewed from a statistical signal processing perspective. Each antenna array 
5 output represents a phase and amplitude coherent sensor which contains a weighted 
sum of the desired signal, the undesired signals, and noise. By applying an 
appropriate complex weight to each sensor signal and then summing the outputs, it is 
possible to cancel undesired interference and enhance the desired signal above the 
noise. The sensor weights are described by an equalizer weight vector w. If the 

10 weight vector is adapted in real time in an optimal manner, it is possible to achieve 
performance that is far superior to fixed beam systems. To implement a fully 
adaptive array, a bank of phase and amplitude coherent digital RF receivers is 
required. Recent technology advances in wide band digital RF converter technology 
make such systems not only feasible, but also economical. An Adaptive spatial 

15 equalizer conceptual diagram is illustrated in Figure 3. 
3.1. Algorithm Cost Function 

In order to adapt the spatial equalizer weight vector, it is necessary to define 
an equalizer signal quality optimization function. Here, we consider two possible 
cost functions, MMSE and maximum SI NR. We will show that for slowly varying 

20 channels, these cost functions are equivalent to within a scalar constant. We will 

also show that for rapidly varying channels, with a fixed equalizer weight applied to 
a finite block of received data, the two cost functions result in very different 
solutions and the maximum SINR criteria is often a more sensible choice. First, the 
slowly varying channel case is considered. Then the case where the channels fade 

25 rapidly with respect to the equalizer adaption rate is discussed. 

Under the slowly varying channel condition, if perfect knowledge were 
available for the instantaneous desired channel, the instantaneous co-channel 
interference channels, the desired signal variane, the interference signal variances, 
and the noise variance, the well known "instantaneous" MMSE spatial equalizer 

30 solution is 

r . min_ ft • . ->) 2 

AW.* = **sr w E [\w x k -s 0 k \-j = o ao R xlt h 0k (7) 
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where R x k is the "instantaneous covariance" matrix for the received signal vector 
defined by 

and R„ k is the instantaneous interference plus noise covariance matrix defined by 

*u,*=£ °h h j t k h L^O Z n I (9) 



The above expressions are derived under the assumption that the desired signal, all 
interfering signals, and the noise are uncorrected. 

Given the same apriori knowledge, the maximum SINR problem is 

- ara 4 "X * g o. J 2 } 
*W - ar 9 w liX p Z M-r (10) 



The solution for W SNR k is easily derived as the generalized eigenvector associated 
with the maximum generalized eigenvalue of the matrix pencil 



l h 0,.* h O..Jf" R u..*l 



which is equivalent to the solution 



W SNR.k " R u!k h 0 .k dD 

0 since 

is obviously rank 1. Using the matrix inversion lemma, we can prove that the 
MMSE solution to Equation (7) is equivalent to the maximum SINR solution to 
Equation (11) to within a scalar constant given by 



MSE, k 



U2) 



So, we conclude that with perfect instantaneous desired and undesired channel 
5 knowledge the MMSE and SINR cost functions lead to nearly identical solutions. 
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In rapidly time varying channel setting, perfect instantaneous knowledge of 
the channels is not available. The channels often change appreciably over the 
adaptation interval of the equalizer algorithm. In such cases, there can be a 
substantial difference between the MMSE and SINK adaption criteria. Consider an 
adaptive spatial equalization algorithm based on an optimal MMSE solution similar 
to Equation (7), but the algorithm employs a sliding window of length N samples for 
each weight vector. Numerous proposed schemes apply a sliding window. The 
sliding window MMSE problem statement may be written as 



where the received signal vector sequence 



(13) 



10 and the transmitted desired signal sequence 

S 0 e k~ I S 0, k" S 0, k-N+l 1 



Again assuming perfect knowledge of the channels, we can show that the solution is 



IN) 



KkK 



(14) 



15 



where 

is the N-sample averaging operator. We can compare this finite block length 
solution to the ideal instantaneous solution to Equation (7) by re-writing Equation 
(14) as 



(15) 



where the covariance error term is 



(16) 



- 19- 



<WO 9B18271A2 I > 



WO 98/18271 



PCT/US97/18196 



and the channel error term is 



* h *4JL h °-~^ ho '* (17 > 



For a time varying channel, as N increases, the error terms in the finite block length 
solution increase. The effect of the covariance error AR k is to add false components 
to the interference subspace in the equalizer solution. The effect of the channel 
error Ah k is additive. The combination of these effects is highly non-linear. Thus, 
it is not clear that applying the MMSE criteria to finite length time blocks with 
rapidly fading channels will result in desirable performance. As a simple example, 
given the zero mean channel model from the previous section, it can be shown that 
lim (N) - 

j^ 00 w MSE.* = 0 (18) 

which is clearly a useless solution. This motivates a search for cost functions which 
are more robust in time varying channels. 

In contrast to the above, applying the SINR criteria to a block of received 
data always results in a useful equalizer solution, even in the limit. Carrying 
through a similar development for the SINR equalizer, we have 



(N) ^ max £||w th 0 , k s o tk "'ho f k-N+i s o,k-N+i)\\l) 



<19) 



The solution for 



w m k 

W SNR ' A 



is the generalized eigenvector associated with the maximum generalized eigenvalue 
of the matrix pencil 

<v kK *v i ») N + o 2 n i 



This solution has a precise meaning: "Given N received samples and one equalizer 
weight vector setting, maximize the ratio of the average received signal energy 
divided by the average noise plus interference energy." The theoretically achievable 
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SINR will decrease as N increases because the interference and desired signal 
covariance matrices will become more well conditioned thus reducing the maximum 
eigenvalue of the matrix pencil. However, this degradation is graceful and is easily 

Equalizer Solution 

(a) 
(b) 
(c) 

analyzed. This will be illustrated in Section 4. 

It is also possible to use convenient combinations of instantaneous and 
average channel information to form variations on the SINR criteria. Three 
variations of the SINR equalizer are presented in Table 3.1. As we will see in the 
next section, it is also sometimes advantageous to replace the averaging operation 

<<*>*>„ 

with the exponentially weighted averaging operator given by 

<(*>*>x=£ **' n l-) n , 0<X<1 (20) 

3.2 Algorithm Description 

The highest potential equalizer performance is obtained with the SINR cost 
function based on the instantaneous desired channel and the instantaneous undesired 
interference subspace. A second solution is possible which is based on a short term 
time average interference subspace description and the instantaneous desired 
channel. The performance for this equalizer is lower than the first since the 
interference matrix includes subspace components from past interference channel 
states. The relative performance decrease is dependent on the angle of arrival spread 
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and the product of the time window length and the Doppler frequency. A third 
equalizer solution involves subspace averaging for both the desired and undesired 
channels with still lower performance. 

From the above, it would seem obvious to choose the instantaneous criteria. 
However, in time varying interference dominated channels, blind equalization 
algorithms which rely on signal property restoral are subject to false signal tracking. 
False tracking occurs when the instantaneous desired channel is in a deep fade and 
the algorithm converges to a stronger interference signal. With a time average 
interference covariance solution, the algorithm becomes less susceptible to false 
tracking as the effective interference subspace averaging window is increased. 
Thus, a trade-off exists between achieving the optimum instantaneous SINR 
performance and robust signal tracking. Jhis motivates us to choose an SINR 
equalizer solution which has the form 

(21) 



where (•) implies an estimated value. We note that the trade-off between 
instantaneous SINR performance and tracking robustness can easily be adjusted by 
varying the forgetting factor A. Increasing X results in a decrease in the rate of false 
tracking occurrences, and a decrease in performance as compared to the ideal 
instantaneous solution. Several techniques are available to detect and recover from 
false tracking events, so signal quality is maintained as long as the occurrence rate is 
low. The proposed equalization strategy is divided into five steps which are 

1. Estimate the present signal sequence by passing the equalizer output 
through a decision device. [Equations (22) and (23)]. 

2. Estimate the present channel state with an optimal estimation filter. 
[Equations (25) and (27)]. 

3. Estimate the received interference plus . noise signal sequence.. Use this 
error signal to update the interference covariance estimate. {Equations (30) and 
(31)]. 

4. Update the spatial equalizer weight vector given the estimated desired 
channel and interference covariance. [Equation (21)]. 
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5. Repeat the process with new desired signal estimates. 
This adaptation procedure is represented by the block diagram shown in Figure 3. 
Each of the algorithm steps are described below. 

Using the present equalizer setting, a short block of received data is 
5 equalized 

^r^X^V'^) (22) 

The equalizer output is then passed through a non-linear detection function d(-) 
which restores some known property of the transmitted portable signal. In the case 
of an FM signal set, the nonlinear function can simply force the signal modulus to a 
constant value, or a more sophisticated combination modulus forcing and phase 
10 filtering function has been proposed. In the case of finite alphabet digital signals, 
the nonlinear function is a symbol waveform detector. The output of the detection 
function is the estimated signal sequence 

5 o,* = d < 5 o,J (23) 

In order to accurately determine the interference siibspace and the optimal 
weight vector, it is critical to have an accurate estimate of the desired channel. The 
15 channel exhibits highly structured time auto-correlation behavior for each antenna 
element. This structured variation should be exploited to improve our ability to 
identify and track the desired channel. The channel is estimated in two steps. First, 
a "measurement" of the reverse channel h O K is obtained by solving 

Here, the channel is assumed time invariant over the data block. The data block 
20 must be a very small fraction of the inverse Doppler frequency to maintain this 

assumption. Assuming the noise plus interference is Gaussian, the solution to 
Equation (24) is the deterministic maximum likelihood estimate 

The channel measurements defined by Equation (25), provide a true channel plus 
noise measurement model since 
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Jtti o.* = V* tJ '* i o.* = V** u * (26) 

where N r = [n lt —n K . N+1 ] and u r is the equivalent channel measurement noise term 

The final channel estimate is a weighted sum of the present and past channel 
measurements 

h* = gf n v n (27) 

n-0 



with the FIR filter vector f given by 

f^rg^E{\\h k - ho J\^ (28) 

5 The solution to Equation (28) is a Wiener filter given by 

where 

is a vector of scalar fading variables representing the time varying complex 
amplitude of one channel vector entry in h K and ^ K = 

is the corresponding time varying measurement noise vector for each entry in u K . 
10 The filter solution to Equation (29) can be based on ensemble average statistics or 
based on worst case sample path behavior. The ensemble average covariance 
function for p K is: 

The noise covariance is determined by the nature of the additive receiver noise and 
the short term cross-correlation and autocorrelation behavior of the interference 
15 signals. 
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The solution for the channel estimation filter is also dependent on two 

■> 

SNR ( — ) 

■? 

°z 

variables, the maximum Doppler shift and the channel measurement . In the interest 
of complexity reduction, the filter design can be based on the maximum expected 
Doppler shift and the minimum expected SNR. Performance is enhanced if the 
Doppler and SNR are estimated and the filter coefficients are adapted either by 
resolving Equation (29) or by choosing the estimation filter from a fixed set of pre- 
computed possibilities. 

To estimate the undesired interference subspace, the undesired signal 
sequence is estimated 

-A^o,* (30) 

The interference covariance estimate is then updated 

<Oa =X *„,*m + <*<* (31) 

The weight vector is now given by Equation (21). The sample index k is 
incremented by N and the next iteration is initiated. 

It should be noted that many variations on this basic algorithm are possible. 
Several implementations involving sample by sample recursions have been 
evaluation. The FIR channel estimation filter can be replaced by an observer based 
filter with various state-space model possibilities. The fundamental contributions are 
the blind SINR adaption criteria, the subspace decomposition and the use of 
correlated desired channel behavior to improve estimation accuracy for the desired 
channel and the interference subspace. 

4. SIMULATIONS 

Simulations have been conducted to test the spatial equalizer approach 
discussed in the previous section. The signals are narrow-band FM with the 
following parameters: 

f 0 = 1GHz (RF center frequency) 
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f D = ± 103Hz (v = 70mph) 

B voicc = 4KHz (voice waveform bandwidth) 

Afvoicc = _±12KHz (peak voice FM excursion) 

fsAT = 6KHz (SAT tone offset) 

5 Af SAT =■ ±2KHz (peak SAT tone FM excursion) 

Here, we present one example simulation scenario in which the base station 
antenna is a six element linear array with inter-element spacing of A/2. There are 
three transmitted mobile signals, each with equal average strength. Each signal 
arrives from four different paths. Each of the four paths for each signal are spread 
10 over a 30° range. Each path for each signal has independent fading generated by 10 

randomly located local scatterers. The central path angles of arrival for the three 
signals are 40°, 20° and -70°. Complex Gaussian band limited noise is added at the 
output of each antenna element. The ratio of the time averaged received signal 
power to average additive noise power is 20dB per element. As the instantaneous 
15 channel fades, the instantaneous signal to additive noise power at each element 

varies above and below 20dB. 

The blind spatial equalization algorithm described in the previous section was 
used to recover estimates for each of the three FM signals. The detection function 
d(-) is the complex sign. The data block length NT is 200/iS with a corresponding 
20 adaptation rate of 5KHz. This is sufficient to support the assumption that the 

channel is constant over each block. The blind algorithm restarts when false signal 
tracking events are detected. False tracking detection can be implemented by several 
methods including monitoring SAT tone quality. 

The instantaneous SINR for the recovered signal waveforms Sq M at the output 
25 of the spatial equalizer for each of the three signal estimates is plotted (solid) along 

with the corresponding optimum theoretical instantaneous SINR performance (dash) 
in Figures 4 through 6. The instantaneous SINR for one omni-directional antenna 
element (dash-dot) is also plotted to provide a reference for SINR improvement. 
Optimum performance is determined by the maximum eigenvalue of the 
30 instantaneous matrix pencil 
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All three simulation runs were performed on the same data set with the same initial 
weight vector setting. The only difference is the signal selected for tracking. 
Conventional CM algorithms were also applied to the simulated data set. None of 
the CM algorithms were found to converge to any of the three signals. 
5 The one second time average estimated signal SINR for signals 1 through 3 is 

18.8dB, 17.2dB and 15.7dB. The average estimated signal SINR loss with respect 
to the optimal instantaneous solution is 0.7dB, 1.9dB, and 3.8dB respectively. This 
is accounted for by the performance loss associated with the subspace forgetting 
factor X. Figure 4 illustrates the theoretical limits for instantaneous SINR 
10 performance sacrifice for signal 3 due to increasing the interference subspace 

forgetting factor X. In Figure 4, the optimal instantaneous SINR solution (solid) is 
plotted along with the maximum eigenvalue of the pencil 

(dash) with X = 0.9999 which approaches infinite averaging. The theoretical 
performance penalty, averaged over one second, between the optimal instantaneous 
15 interference covariance solution and the infinite average interference covariance 

solution is 4.6dB for the signal 3 example. These simulation results suggest that the 
proposed algorithm achieves optimal theoretical performance when the theoretical 
penalty due to interference subspace averaging is taken into account. 

CONCLUSIONS 

20 A new approach for blind adaptive spatial equalization in rapidly fading 

channels has been developed. The algorithm performs well in an extreme narrow 
band FM communication environment with high Doppler spread, high angle spread, 
and strong multiple co-channel interferers. We believe that these are the first 
reported results for a blind spatial equalizer which performs well on AMPS type 

25 signals in a high velocity portable multipath setting. The approach involves two 
fundamental results: 

1. For certain communication signals in rapidly fading channel settings, a 
modified version of the SINR adaptive equalizer cost function is superior to 
conventional squared error cost functions. 
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2. Optimal estimation techniques can be applied to obtain accurate channel 
estimates in the presence of noise and interference. High quality channel estimates 
are crucial to obtaining an accurate estimate of the interference subspace used to 
determine the SINR equalizer solution. 

Extensions of our approach have been developed for multiple antenna delay- 
spread digital FIR communication channels. These time varying space-time 
equalizer results will be reported in a future paper. 
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WHAT IS CLAIMED IS: 

1. A base station for communicating over a radio network with a remote station, 
comprising: 

a plurality of antennas and receivers to receive inbound information signals 
from the remote station; 

a signal processor coupled to said receivers and configured to compute a 
plurality of weights for the purpose of weighting each of the receiver information 
signals and combining said weighted signals together in a manner that improves the 
signal quality at the output of said signal processor for both geographically 
stationary and rapidly moving remote stations. 

2. The base station of claim 1, wherein said signal processor includes: 
means for statistically characterizing a receive channel vector, said receive 

channel vector being representative of a receive communications channel over which 
signal energy is transmitted from said at least one remote communications station to 
said central communications station; and 

means for generating said plurality of weights based on results of statistical 
characterization of said receive channel vector. 

3. The base station of claim 2, wherein said signal processor includes: 
one or more time-varying channel estimation filters for estimating said 

receive channel vector. 

4. The base station of claim 2, wherein said means for statistically 
characterizing further include means for compiling an estimated receive channel 
covariance matrix. 

5. The base station of claim 4, wherein said means for generating said beam 
pattern weight vector includes means for averaging said estimated receive channel 
covariance matrix in order to form an estimated transmit channel covariance matrix. 
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6. The base station of claim 5, further including means for using said estimated 
transmit channel covariance matrix to form objective and constraint functions, said 
beam pattern weight vector being obtained by solving said objective and constraint 
functions. 

7. A method of communicating over a radio network with a remote station 
using a base station having a plurality of antennas and receivers and a signal 
processor, said method comprising the steps of: 

receiving inbound information signals from the remote station via said 
plurality of antennas and receivers; 

computing a plurality of weights in said signal processor for the purpose of 
weighting each of the receiver information signals and combining said weighted 
signals together in a manner that improves the signal quality at the output of said 
signal processor for both geographically stationary and rapidly moving remote 
stations. 

8. The method of claim 7, wherein said computing step includes the steps of: 
statistically characterizing a receive channel vector, said receive channel 

vector being representative of a receive communications channel over which signal 
energy is transmitted from said at least one remote communications station to said 
central communications station; and 

generating said plurality of weights based on results of statistical 
characterization of said receive channel vector. 

9. The method of claim 8, wherein said method further includes the step of: 
estimating said receive channel vector. 

10. The method of claim 8, wherein said statistically characterizing step includes 
the step of: 

compiling an estimated receive channel covariance matrix. 
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11. The method of claim 10, wherein said generating said beam pattern weight 
vector step includes the step of: 

averaging said estimated receive channel covariance matrix in order to form 
an estimated transmit channel covariance matrix. 

12. The method of claim 1 1 , further comprising the step of: 

using said estimated transmit channel covariance matrix to form objective and 
constraint functions, said beam pattern weight vector being obtained by solving said 
objective and constraint functions. 
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WIRELESS COMMUNICATION NETWORK USING TIME. VARYING 
VECTOR CHANNEL EQUALIZATION FOR ADAPTIVE SPATIAL 

EQUALIZATION 

RELATED APPLICATIONS 
5 This application claims priority to the following applications and incorporates 

these applications by reference: 

U.S. Ser. No. 08/394,652 filed on February 22, 1995; 
U.S. Ser, No. 08/491,044, filed on June 16, 1995; and 
U.S. Prov. App. No. 60/005,647 filed on October 19, 1995. 

10 FIELD 

The present invention relates to a wireless communication network that uses 
time- varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
quality of communication. 

15 BACKGROUND 

Within wireless mobile communication systems, procedures for compensating 
for multipath effects are desirable. Unlike "line-of-sight" radio links, a number of 
signal transmission paths typically comprise each wireless communication channel, 
hence the term "multipath." An increase in primary path communication energy is 

20 desired with an attendant reduction in the interference energy radiated to mobile 
users over other non-primary paths. Often this increase in primary path 
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communication energy is achieved through generation of spatially selective, directive 
transmission beam patterns. 

Some systems employ directive antennas at base station sites to increase the 
signal level received from and transmitted to each mobile user relative. Other 
5 systems use an antenna array that is used to form beams to increase the signal level 
received from and transmitted to each mobile user. In the systems with an antenna 
array, often called a phased array, weights are used to account for various angles 
that the beam will provide. Problems occur when severe multipath signals are 
created by a plurality of obstructions such as buildings and mountains. When severe 
10 multipath occurs, many systems suffer severe quality degradation. Known systems 
cannot sufficiently account for severe multipath and the result is reduced quality 
communication. 

Figure 1 shows an illustrative representation of a wireless "multipath" 
communication channel between a base station 2 and a remote mobile user 4. The 

15 various signal transmission paths within each such multipath communication channel 
arise from reflection of the transmitted signal by dominant reflection surfaces 6, and 
by minor reflection surfaces 12, between the base station 2 and remote mobile user 
4. Accordingly, techniques for improving signal reception in line-of-sight radio 
systems are often not directly applicable to multipath signal environments. For 

20 example, in line-of-sight system the "gain" of an antenna typically varies inversely 
to antenna beam width. However, if a given antenna beam width is made less than 
the angular region encompassing the various signal paths comprising a multipath 
communication channel, further reduction in the beam width may reduce the energy 
radiated along some of the angular paths. In turn, this may actually decrease the 

25 effective time average gain of the antenna. 

These systems that perform blind transmission fail to employ a sufficient 
method for determining the received angle of the mobile station signal. As a result, 
known transmission beamforming techniques often use a trial and error approach 
that works barely satisfactorily, and does not always provide high quality 

30 communication. 
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RELATED ART 

Within wireless mobile communication systems, three techniques have been 
developed for improving communication link performance using directive transmit 
antennas: (i) selection of a particular fixed beam from an available set of fixed 
5 beams, (ii) adaptive beam forming based on receive signal angle estimates, (iii) 

adaptive transmission based on feedback provided by the remote mobile user, and 
(iv) adaptive transmit beam forming based upon the instantaneous receive beam 
pattern. Each of these techniques is described briefly below. 

In the first technique, one of several fixed base station antenna beam patterns 

10 is selected to provide a fixed beam steered in a particular direction. The fixed 
antenna beams are often of equal beam width, and are often uniformly offset in 
boresight angle so as to encompass all desired transmission angles. The antenna 
beam selected for transmission typically corresponds to the beam pattern through 
which the largest signal is received. The fixed beam approach offers the advantage 

15 of simple implementation, but provides no mechanism for reducing the signal 

interference power radiated to remote mobile users within the transmission beam of 
the base station. This arises because of the inability of the traditional fixed beam 
approach to sense the interference power delivered to undesired users. 

The second approach involves "adapting" the beam pattern produced by a 

20 base station phase array in response to changing multipath conditions. In such 
beamforming antenna arrays, or "beamformers", the antenna beam pattern is 
generated so as to maximize signal energy transmitted to ("transmit beamforming"), 
and received from ("receive beamforming"), an intended recipient mobile user. 
While the process of transmit beamforming to a fixed location over a 

25 line-of-sight radio channel may be performed with relative ease, the task of 

transmitting to a mobile user over a time-varying multipath communication channel 
is typically considerably more difficult. One adaptive transmit beamforming 
approach contemplates determining each angle of departure (AOD) at which energy 
is to be transmitted from the base station antenna array to a given remote mobile 

30 user. Each AOD corresponds to one of the signal paths of the multipath channel, 
and is determined by estimating each angle of arrival (AOA) at the base station of 
signal energy from the given user. A transmit beam pattern is then adaptively 
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formed so as to maximize the radiation projected along each desired AOD (i.e, the 
AOD spectrum), while minimizing the radiation projected at all other angles. 
Several well known algorithms (e.g., MUSIC, ESPRIT, and WSF) may be used to 
estimate an AOA spectrum corresponding to a desired AOD spectrum. 
5 Unfortunately, obtaining accurate estimates of the AOA spectrum for 

communications channels comprised of numerous multipath constituents has proven 
problematic. Resolving the AOA spectrum for multiple co-channel mobile units is 
further complicated if the average signal energy received at the base station from any 
of the mobile units is significantly less than the energy received from other mobile 

10 units. This is due to the fact that the components of the base station array response 
vector contributed by the lower-energy incident signals are comparatively small, thus 
making it difficult to ascertain the AOA spectrum corresponding to those mobile 
units. Moreover, near field obstructions proximate base station antenna arrays tend 
to corrupt the array calibration process, thereby decreasing the accuracy of the 

15 estimated AOA spectrum. 

In the third technique mentioned above, feedback information is received at 
the base station from both the desired mobile user, and from mobile users to which 
it is desired to minimize transmission power. This feedback permits the base station 
to "learn" the "optimum" transmit beam pattern, i.e., the beam pattern which 

20 maximizes transmission to the desired mobile user and minimizes transmission to all 

other users. One disadvantage of the feedback approach is that the mobile radio 
needs to be significantly more complex than would otherwise be required. 
Moreover, the information carrying capacity of each radio channel is reduced as a 
consequence of the bandwidth allocated for transmission of antenna training signals 

25 and mobile user feedback information. The resultant capacity reduction may be 

significant when the remote mobile users move at a high average velocity, as is the 
case in most cellular telephone systems. 

The fourth conventional technique for improving communication link 
performance involves use of an optimum receive beam pattern as the preferred 

30 transmission beam pattern. After calibrating for differences between the antenna 
array and electronics used in the transmitter and receiver, it is assumed that the 
instantaneous estimate of the nature of the receive channel is equivalent to that of the 
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transmit channel. Unfortunately, multipath propagation and other transient channel 
phenomenon can substantially eliminate any significant equivalence between 
frequency-duplexed transmit and receive channels, or between time-division 
duplexed transmit and receive channels separated by a significant time interval. As 
5 a consequence, communication link performance fails to be improved. 

SUMMARY 

The present invention relates to a wireless communication network that uses 
time-varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
10 quality of communication. 

Accordingly, it is an object of the invention to provide an adaptive receive 
and transmit beamforming technique that improves communication quality by 
utilizing the received signal energy from remote users. 

It is another object of the invention to provide an adaptive receive and 
15 transmit beamforming technique which accounts for the presence of multipath 

fading inherent in the communication channel. 

It is yet another object of the invention that the beamforming technique be 
independent of antenna array geometry, array calibration, or of explicit feedback 
control signals from remote users. 
20 It is another object of the invention to provide adaptive receive and transmit 

beam forming that improves signal quality received by a desired user and while 
simultaneously reducing interference energy received by other undesired users so as 
to, within a cellular communication network, improve communication traffic 
capacity, and/or to increase base station coverage area, and/or to improve call 
25 quality. 

An adaptive reception and transmission technique according to one 
embodiment of the invention offers the advantages of adaptive reception 
transmission using feedback without the associated mobile radio complexity increase 
and information capacity penalty. The technique has been developed to exploit 
30 structured variation which occurs in the multipath fading present in the wireless 

antenna array channel. Thus, multipath propagation effects are explicitly accounted 
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for in the problem approach. The technique is blind in that the antenna beam is 
formed in the absence of explicit knowledge of the array geometry, and without the 
necessity of array calibration or mobile feedback. The basic approach is to estimate 
the optimum receive and transmit antenna beam pattern based on certain statistical 
5 properties of the received antenna array signals. The optimum receive and transmit 

beam pattern is found by solving an optimization equation. The adaptive transceiver 
system is suitable for use in conjunction with either a diplexed transmit/receive 
antenna array, or with separate transmit and receive arrays. 

Advantages of the invention include reduced interference and improved 
10 communication quality. 

BRIEF DESCRIPTION OF THE FIGURES 
Additional advantages of the invention will become apparent upon reading 
the following detailed description and upon reference to the drawings, in which: 

Figure 1 is an illustrative representation of a multipath propagation channel 
15 between a base station and a remote user station; 

Figure 2A is a block diagram of the physical implementation of a 
beamforming network configured to perform adaptive beam forming in accordance 
with the present invention; 

Figure 2B is a block diagram of the physical implementation of a 
20 beamforming network in which a diplexor is employed to allow the antenna array to 
be used for both transmit and receive operation; 

Figure 3 is a functional block diagrammatic representation of a beamforming 
network of the present invention; 

Figure 4 is a graph depicting estimated signal SINR versus adaptation time 
25 for signal 1 ; 

Figure 5 is a graph depicting estimated signal SINR versus adaptation time 
for signal 2; 

Figure 6 is a graph depicting estimated signal SINR versus adaptation time 
for signal 3; and 

30 Figure 7 is a graph showing theoretical performance bounds for variable X. 
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DETAILED DESCRIPTION 
The present invention relates to a wireless communication network that uses 
time-varying vector channel equalization for adaptive spatial equalization. In 
particular, the invention is used in a wireless communication network to improve the 
5 quality of communication. 

The exemplary embodiments are described herein with reference to specific 
configurations and protocols. Those skilled in the art will appreciate that various 
changes and modifications can be made to the exemplary embodiments while 
remaining within the scope of the present invention. 

10 The problem we consider is blind adaptive spatial equalization with a 

multiple antenna base station array for the purpose of increasing reverse link spectral 
efficiency (Erlangs/Hz/base station) or range by canceling interference and 
enhancing the desired signal over noise. In blind equalization, the equalizer weight 
vector adaptation relies on the structure of the transmitted signal rather than a known 

15 training signal. While these techniques may be applied at the portable unit, this 

present embodiments focus on spatial equalization at the base station. 

The communication channel formed between a portable wireless user and a 
wireless base station varies with time. When the portable moves with high velocity, 
the time variation in the channel can become severe. Over the past decade, a large 

20 body of work has been devoted to blind spatial equalization in time invariant 
channels. However, less work has been reported for large Doppler-spread 
conditions. For certain communication signals, including narrow-band FM, the 
known existing equalizer algorithms fail to converge in rapidly fading environments. 
This occurs because the channel changes significantly over the time interval required 

25 to determine the next equalizer weight vector update. 

In most conventional adaptive algorithms, the signal quality criteria is 
minimum mean squared error (MMSE). We will show that this criteria may not be 
the best choice in rapidly fading channels. We will also show that the signal to 
interference plus noise ratio (SINR) criteria can be more robust in rapidly fading 

30 channels. With this information and motivation, an algorithm is developed which 
adaptively optimizes an interesting variation of the conventional SINR criteria. 
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The exemplary equalizer decomposes the received vector signal into a desired ' 
signal subspace and an interference plus noise subspace. In order to accurately 
estimate these two subspaces, we exploit the inherent time correlation structure that 
exists in the fading channel to obtain instantaneous blind estimates of the desired 
channel. Supplement 1 is provided for additional information and includes 
additional background and introduction material. 



I. O vervi ew of Beamforming Network 

Turning now to Figure 2 A, a block diagram is shown of the physical 
organization of a beamforming network 50 configured to perform adaptive beam 

10 forming in accordance with the present invention. In an exemplary embodiment the 
beamforming network 50 is disposed within a base station of a cellular 
communications network, in which is included a transceiver comprised of a radio 
frequency (RF) transmitter 52 and an RF receiver 54. 

In the embodiment of Figure 2 A, a base station antenna array 56 serves to 

15 produce independent transmit and receive antenna beams for facilitating 

communication with one or more mobile units (not shown). The term "receive 
channel vector" is employed to indicate that each antenna element within the base 
station antenna array 56 will form a propagation channel to a given remote user. 
The composite array channel may be represented as a vector having elements 

20 corresponding to each individual antenna channel. As is described herein, statistical 
characterization of the receive channel vector provides information which may be 
used by the base station to determine an "optimal" transmit beam pattern, i.e., a 
transmit beam pattern which maximizes the average signal-to-interference power 
delivered to a given mobile user. This obviates the need for the mobile unit to 

25 provide feedback information to the base station relating to propagation 

characteristics of the transmit channel. This in turn simplifies implementation of the 
mobile unit, and preserves channel resources for data transmission by eliminating the 
need for "mobile unit feedback" relating channel characteristics to the base station. 
As is indicated by Figure 2A, a diplexer 58 can be employed to allow the 

30 antenna array 56 to be used for both transmit and receive operation by isolating the 

RF receiver 54 from the RF transmitter 52. A receive channel beamformer 60 
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cooperates with the RF receiver 54 to adaptively optimize the receive antenna beam 
in order to improve received signal quality. Similarly, a transmit channel 
beamformer 64 cooperates with the RF transmitter 52 to adapt the transmit antenna 
beam to optimize some characteristic of transmission quality. In an exemplary 
5 embodiment the transit channel beamformer 64 and receive channel beamformer 60 
are each implemented as a special purpose digital signal processor (DSP). 

In another embodiment of the invention, distinct antenna arrays are used for 
signal reception and transmission as illustrated in Figure 2B. In the embodiment of 
Figure 2B, a diplexer is not required since a dedicated receive antenna array (not 

10 shown) is connected to the RF receiver 54, and a dedicated transmit antenna array 

(not shown) is connected to the RF transmitter 52. The receive and transmit antenna 
arrays are designed to provide identical radiation characteristics when operated at the 
receive and transmit frequencies, respectively. Accordingly, in many instances the 
physical geometries of the transmit and receive antenna arrays are simply physically 

15 scaled to account for the fractional difference in the receive and transmit RF 

wavelengths. The embodiment of Figure 2B substantially eliminates the potential 
introduction of error arising from use of a single antenna array and diplexer. 

Turning now to Figure 3, a functional block diagrammatic representation is 
provided of a beamforming network of the present invention. In Figure 3, solid 

20 lines are used to represent functional elements and dashed lines are employed to 
identify the physical components of Figure 2. The RF receiver 54 is functionally 
comprised of a low-noise amplifier (LNA) network 110, and a complex 
downconverter 1 12. The complex downconverter 112 is disposed to frequency 
downconvert the received RF signal energy after processing by the LNA network 

25 1 10. The downconverted signal energy is then digitally sampled and provided to a 

receive channel weighting module 1 16 of the receive channel beamformer 60. The 
weights applied by the receive channel beamformer 60 to each of the M 
downconverted antenna element outputs x R ni (k), m =1 to M, of the complex 
frequency downconverter 1 12 are determined by a receiver weight vector adaptation 

30 module 118. In the exemplary embodiment the receiver weight vector adaptation 
module 118 determines a receive channel weight vector, w R , which maximizes the 
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signal quality received over the desired inbound frequency channel. The weight 
vector calculation is performed according to that described in Supplement 1. 

In the embodiment of Figure 3, a vector channel covariance estimator 140 
within the transmit beamformer 64 operates to produce a statistical characterization 
5 of a receive channel vector using: (i) the outputs x Rm (k), m = 1 to M, of the 

complex frequency downconverter 1 12, and (ii) an estimate of the desired signal 
§ r (k) generated in the receive channel beamformer 60. For present purposes the 
receive channel vector may be viewed as being representative of the multipath 
communications channel from a mobile user (not shown in Figure 3) to the antenna 

10 array 56. In an exemplary embodiment the statistical characterization carried out 
within the covariance estimator 140 yields an estimated receive channel covariance 
matrix used during the transmit beamforming process. Throughout the following 
description, scalar quantities are represented using lower case characters of standard 
type, vector quantities are represented using lower case characters of bold type, and 

15 matrix quantities are represented using upper case characters of bold type. 

Within the transmit channel beamformer 64, an optimal transmit beam 
pattern weight vector is generated by a transmit channel weight vector adaptation 
module 150 based on the results of the statistical characterization of the receive 
channel vector. This weight vector, w T (t), optimizes a particular aspect (e.g., 

20 average desired signal to undesired interference ratio) of the signal energy within the 
transmission range of the base station array 56. In the exemplary embodiment, the 
optimal transmit beam pattern weight vector is generated using the estimated desired 
receive channel covariance matrix, R^k), and undesired interference covariance 
matrix, R^k), both of which are compiled within the covariance estimator 140. 

25 As is indicated by Figure 3, the signal information (S T ) to be transmitted to 

the desired mobile radio unit is used to modulate a digital baseband carrier within a 
modulator 154. The modulated signal is then applied to a transmit channel 
weighting module 158 disposed to weight, on the basis of the optimized transmit 
pattern weight vector, the input signals corresponding to each element of the 

30 antenna array 56. The weighted set of input signals produced by the weighting 
module 158 are then upcon verted in frequency within a complex frequency 
upconverter 160 of the RF transmitter 52. The resultant frequency-upconverted 

- 10- 



^nonin- <WO__ 981 8271 A3 I A- 



WO 98/18271 



PCT/US97/18196 



signals are then amplified by a bank of power amplifiers 164, and provided to the 
antenna array 56 for transmission via diplexer 58. 

In the exemplary embodiment an improved estimate of the received signal is 
obtained through utilization of a demodulator 170 and remodulator 172. The 
5 received signal is demodulated within demodulator 170 in order to recover the 
essential characteristics of the modulating signal. In the case of an analog FM 
signal, this involves recovery of the FM waveform. In the case of a digitally 
modulated signal (BPSK, FSK, QPSK, etc.), the demodulator 170 forms hard 
decisions as to the value of each digital symbol. The demodulated signal is then 

10 processed based upon some predefined characteristic of the signal and modulation 
pattern. For example, a demodulated analog FM signal could be lowpass filtered 
based upon a known signal information bandwidth as a means of obtaining an 
improved post-demodulation signal estimate. In the case of digital modulation, error 
correction could be implemented in order to remove bit errors, thus improving 

15 estimated signal quality. In addition, training signals (i.e., pilot tones, SAT tones, 

etc.) may optionally be employed in lieu of, or in conjunction with, the 
aforementioned "blind" techniques. 

Again referring to Figure 3, the processed demodulated signal is then used by 
demodulator 172 to remodulate an RF carrier, thereby producing an improved 

20 modulated signal estimate. The improved signal estimate is then used by the 

receiver weight vector adaptation block 1 18 and the covariance estimate 140. Other 
techniques, which do not rely upon such a demodulation/remodulation procedure, 
can be devised for obtaining a sufficiently accurate received signal estimate. Figure 
3 simply illustrates a particular exemplary embodiment incorporating a demodulator 

25 and remodulator. 

In the present embodiment, the demodulator 170 and remodulator 172 or the 
receive channel beamformer 60 are operative to produce a received signal estimate 
s r . The quantity S r is then employed by the covariance estimator 140 to estimate the 
covariance matrix of the receive channel. The receive channel beamformer 60 could 

30 of course be replaced by numerous alternative structures including, for example, 

multi-antenna sequence estimators, maximum likelihood receiver structures, multiple 
element decision feedback equalizers, or similar alternative described in Supplement 
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1. Any of these alternative structures could also be used to provide the quantity s r 
for use in estimating the received channel covariance statistics. 

IL Blind Beamforming Using Multipart! Signals 
As mentioned above, an initial step within the present invention involves 
5 characterizing a receive channel vector representative of the multipath 

communications channel from a mobile user to a base station antenna array 56. In 
one embodiment, this is accomplished as described in Supplement 1 . This 
embodiment employs a five step approach: 

1) estimate the present signal sequence by passing the equalizer output 
10 through a decision device; 

2) estimate the present channel state with an optimal estimation filter; 

3) estimate the received interference plus noise signal sequence and use this 
error signal to update the interference covariance estimate; 

4) update the spatial equalizer weight vector given the estimated desired 
15 channel and interference covariance; and 

5) repeat the process with new desired signal estimates. 

III. Alternate Embodiments 

Although the embodiments of the invention presented herein have been 
described in terms of optimization from a base station perspective, the teachings of 

20 the invention are equally applicable to optimum formation of the mobile unit antenna 
pattern. That is, the mobile unit antenna pattern is formed such that a desired level 
of power is transmitted to the one or more base stations with which the mobile unit 
is in communication, and so that transmission power is minimized to the other base 
stations within the communication network. Again, an initial step in this antenna 

25 pattern optimization procedure involves statistically characterizing the mobile unit 
receive channel vector. The results of this statistical characterization are then used 
to determine the beam pattern weight vector which maximizes a predetermined 
quality parameter (e.g., signal to noise ratio) of the signal received by the intended 
base station(s), while minimizing signal transmission to other base stations. At the 
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mobile station, the transmit antenna beam pattern is then generated in accordance 
with the transmit beam pattern weight vector. 

Accordingly, in the following claims the term "central communications 
station" is intended to refer to either the base station or mobile unit configured to 
generate an optimized antenna beam pattern in accordance with the invention. When 
a base station comprises the central communications station, the mobile unit(s) in 
communication therewith are referred to as "remote communications station(s)". 
Conversely, when a mobile unit comprises the central communications station, the 
receiving base station(s) are identified as remote communications station(s). 

While the present invention has been described with reference to a few 
specific embodiments, the description is illustrative of the invention and is not to be 
construed as limiting the invention. Various modifications may occur to those 
skilled in the art while remaining within the true spirit and scope of the invention as 
defined by the appended claims. 

SUPPLEMENT 1 

New blind adaptive spatial equalization techniques are proposed for severe 
multipath fading with co-channel interference. In rapidly fading environments, 
equalizer performance suffers with existing algorithms because the desired and 
interference channels change significantly over the update time interval. It is shown 
that for fast fading channels, conventional squared error equalizer cost functions can 
yield poor results. A new algorithm is proposed based on a variation of the SINK 
cost function coupled with optimal time varying channel estimation to improve 
signal tracking. In simulations, the algorithm achieves theoretically optimum time 
varying equalizer performance in extreme interference channels with narrow band 
FM signals and high velocity portables. These results are contrary to the widespread 
belief that high performance blind adaptive array processing is not possible for 
AMPS with severe combinations of Doppler spread, angle spread, and interference 
power. 
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1. INTRODUCTION 

The problem we consider is blind adaptive spatial equalization with a 
multiple antenna base station array for the purpose of increasing reverse link spectral 
efficiency (Erlangs/Hz/base station) by canceling interference and enhancing the 
desired signal over noise. In blind equalization, the equalizer weight vector 
adaptation relies on the structure of the transmitted signal rather than a known 
training signal. While these techniques may be applied at the portable unit, this 
paper focuses on spatial equalization at the base station receiver. 

The communication channel formed between a portable wireless user and a 
base station varies with time. When the portable moves with high velocity, the time 
variation in the channel can become severe. Over the past decade, a large body of 
work has been devoted to blind spatial equalization in time invariant channels. 
However, less work has been reported for large Doppler-spread conditions. For 
certain communication signals, including narrow-band FM, the known existing 
equalizer algorithms fail to converge in rapidly fading environments. This occurs 
because the channel changes significantly over the time interval required to 
determine the next equalizer weight vector update. 

In most conventional adaptive algorithms, the signal quality criteria is 
minimum mean squared error (MMSE). We will show that this criteria may not be 
the best choice in rapidly fading channels. We will also show that the signal to 
interference plus noise ratio (SINR) criteria can be more robust in rapidly fading 
channels. With this motivation, an algorithm is developed which adaptively 
optimizes an interesting variation of the conventional SINR criteria. The proposed 
equalizer decomposes the received vector signal into a desired signal subspace and 
an interference plus noise subspace. In order to accurately estimate these two 
subspaces, we exploit the inherent time correlation structure that exists in the fading 
channel to obtain instantaneous blind estimates of the desired channel. 

The multiple antenna wireless channel model is briefly reviewed in Section 2. 
The spatial equalization problem is developed and a solution is proposed in Section 
3. Simulation results demonstrating the performance of the adaptive spatial 
equalizer with multiple co-channel narrow-band analog FM signals in a high 
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Doppler spread, high angle spread fading environment are presented in Section 4. 
Conclusions are offered in Section 5. 

2. CHANNEL MODELS 
The wireless communication channel is characterized by scattering and 
5 reflecting bodies which cause severe multipath. This is illustrated in Figure 1. The 
channel is also time- varying when the receiver is in motion. The time- variation is 
characterized by rapid amplitude fluctuations and therefore is called a fading 
channel. These short-term complex amplitude variations are often well modeled by 
a multiplicative Gaussian random process. The channel also exhibits long-term 
10 variations which are due to shadowing by large physical features. This variation is 
much less structured and is sometimes well modeled as a log-normal multiplicative 
random variable. 

Let s Q (t) be the complex analytic baseband representation of the transmitted 
signal. If the remote transmitter has a single antenna and the wireless base station 
15 receiver has multiple antennas, then the direction of arrival of the 1 th path form the 

receiver (given by 6 0 1) is of importance. The baseband vector output of the base 
station receiver dues to the 1 th propagation path is 

x 2 (t)=aiQ 01 )fS 01 (t)s 0 (t) (1) 

where a(6 0 ^ denotes the base station array response to the path arriving at the angle 
8 0i . The composite complex-valued propagation path strength p 0il (t) is given by 

20 where a 0 x (i) is the complex fast fading parameter, T, is the shadow fading loss 

factor, d 0 is the distance from the mobile to the base station and r| is the path power- 
loss exponent. In Equation (1) an inherent assumption is that each propagation path 
amplitude is frequency invariant over the bandwidth of the signal and that the 
relative time delay between different arriving paths is small compared to the inverse 

25 signal bandwidth. Below, it is further assumed that the array response a(6 01 ) is 
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invariant over the signal bandwidth. This set of conditions is commonly referred to 
as the narrow band assumption. 

The received time-domain signal vector is the sum of all dominant reflector 
paths and the additive receiver noise n(t) 

x{t)=^x i (t)+/i{t)-^a(6 0 J {t)s 0 (t)+ii(t) (3) 



Clearly from the above model, the time-varying response of the continuous-time 
desired portable channel is 



The channel given in Equation (3) is sampled at a rate of T seconds. The 
output samples can be rewritten as, 

X k= h o,k s o,k + n k (5) 

with the discrete time channel response given by 

10 In the presence of multiple user co-channel interference, we have the received signal 
vector 

When portables move at high velocities, giving rise to large Doppler 
frequency shifts, the channel vector complex M-space direction and magnitude 
become rapidly varying functions of time. The propagation environment determines 
15 the exact nature of the vector variation. The central path angle and relative spread 

in path angles (angle spread) for each channel determines the Eigenvectors and 
Eigenvalues of the time averaged vector subspace spanned by the channel. The 
Doppler spread determines the rate of time variation within the channel subspace. 
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3. SPATIAL CHANNEL EQUALIZATION 



The techniques proposed in this paper do not consider conventional antenna 
concepts such as beam width, gain and side lobe rejection. Instead, the antenna 
array is viewed from a statistical signal processing perspective. Each antenna array 
output represents a phase and amplitude coherent sensor which contains a weighted 
sum of the desired signal, the undesired signals, and noise. By applying an 
appropriate complex weight to each sensor signal and then summing the outputs, it is 
possible to cancel undesired interference and enhance the desired signal above the 
noise. The sensor weights are described by an equalizer weight vector w. If the 
weight vector is adapted in real time in an optimal manner, it is possible to achieve 
performance that is far superior to fixed beam systems. To implement a fully 
adaptive array, a bank of phase and amplitude coherent digital RF receivers is 
required. Recent technology advances in wide band digital RF converter technology 
make such systems not only feasible, but also economical. An Adaptive spatial 
equalizer conceptual diagram is illustrated in Figure 3. 
3.1. Algorithm Cost Function 

In order to adapt the spatial equalizer weight vector, it is necessary to define 
an equalizer signal quality optimization function. Here, we consider two possible 
cost functions, MMSE and maximum S1NR. We will show that for slowly varying 
channels, these cost functions are equivalent to within a scalar constant. We will 
also show that for rapidly varying channels, with a fixed equalizer weight applied to 
a finite block of received data, the two cost functions result in very different 
solutions and the maximum SINR criteria is often a more sensible choice. First, the 
slowly varying channel case is considered. Then the case where the channels fade 
rapidly with respect to the equalizer adaption rate is discussed. 

Under the slowly varying channel condition, if perfect knowledge were 
available for the instantaneous desired channel, the instantaneous co-channel 
interference channels, the desired signal variane, the interference signal variances, 
and the noise variance, the well known "instantaneous" MMSE spatial equalizer 
solution is 




<7) 
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where R x k is the "instantaneous covariance" matrix for the received signal vector 
defined by 

R X ,S E K**}"°5A>.»Ao. *♦*„., (8) 

and k is the instantaneous interference plus noise covariance matrix defined by 

R u, *= £ a lj h 3. * h /. * + oil ( 9 ) 



The above expressions are derived under the assumption that the desired signal, all 
interfering signals, and the noise are uncorrected. 

Given the same apriori knowledge, the maximum SINR problem is 

- ara max , fl 1 w ' A °.* 5 o,«l1 , 



The solution for W SNR k is easily derived as the generalized eigenvector associated 
with the maximum generalized eigenvalue of the matrix pencil 

which is equivalent to the solution 

10 since 

h 0,k h 0.k 

is obviously rank 1 . Using the matrix inversion lemma, we can prove that the 
MMSE solution to Equation (7) is equivalent to the maximum SINR solution to 
Equation (11) to within a scalar constant given by 



w 

MSE,k 



(12) 



So, we conclude that with perfect instantaneous desired and undesired channel 
15 knowledge the MMSE and SINR cost functions lead to nearly identical solutions. 
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In rapidly time varying channel setting, perfect instantaneous knowledge of 
the channels is not available. The channels often change appreciably over the 
adaptation interval of the equalizer algorithm. In such cases, there can be a 
substantial difference between the MMSE and SINR adaption criteria. Consider an 
adaptive spatial equalization algorithm based on an optimal MMSE solution similar 
to Equation (7), but the algorithm employs a sliding window of length N samples for 
each weight vector. Numerous proposed schemes apply a sliding window. The 
sliding window MMSE problem statement may be written as 

- «g {|| WX k -s Q J\l} (13) 

where the received signal vector sequence 
and the transmitted desired signal sequence 

S 0,k = l S 0,k"' S 0,k-N+l} 

Again assuming perfect knowledge of the channels, we can show that the solution is 



where 



is the N-sample averaging operator. We can compare this finite block length 
solution to the ideal instantaneous solution to Equation (7) by re-writing Equation 
(14) as 

where the covariance error term is 
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and the channel error term is 



(17) 



For a time varying channel, as N increases, the error terms in the finite block length 
solution increase. The effect of the covariance error AR k is to add false components 
to the interference subspace in the equalizer solution. The effect of the channel 
error Ah k is additive. The combination of these effects is highly non-linear. Thus, 
it is not clear that applying the MMSE criteria to finite length time blocks with 
rapidly fading channels will result in desirable performance. As a simple example, 
given the zero mean channel model from the previous section, it can be shown that 



lim 



AN) 



N^ 00 W MS£.k - u (18) 



which is clearly a useless solution. This motivates a search for cost functions which 
10 are more robust in time varying channels. 

In contrast to the above, applying the SINR criteria to a block of received 
data always results in a useful equalizer solution, even in the limit. Carrying 
through a similar development for the SINR equalizer, we have 

tw> A „ niax E V W t h o,k s o,k" h o,k-N+i s o t k~N+i\\\l\ 

SNR.k ~ w —, 7— * ~ r 



The solution for 



W SNR 9 A 



15 is the generalized eigenvector associated with the maximum generalized eigenvalue 

of the matrix pencil 

ko, k K *> w ; i a 2 sj ih^ k hl *> N +a;i 



This solution has a precise meaning: "Given N received samples and one equalizer 
weight vector setting, maximize the ratio of the average received signal energy 
divided by the average noise plus interference energy." The theoretically achievable 
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SINR will decrease as N increases because the interference and desired signal 
covariance matrices will become more well conditioned thus reducing the maximum 
eigenvalue of the matrix pencil. However, this degradation is graceful and is easily 

Equalizer Solution 

(a) 
(b) 
(c) 

analyzed. This will be illustrated in Section 4. 

It is also possible to use convenient combinations of instantaneous and 
average channel information to form variations on the SINR criteria. Three 
variations of the SINR equalizer are presented in Table 3. 1. As we will see in the 
next section, it is also sometimes advantageous to replace the averaging operation 

<<•>*>„ 

with the exponentially weighted averaging operator given by 

£ X *~ n <->„ r 0<X<1 (20) 

n--oo 

3.2 Algorithm Description 

The highest potential equalizer performance is obtained with the SINR cost 
function based on the instantaneous desired channel and the instantaneous undesired 
interference subspace. A second solution is possible which is based on a short term 
time average interference subspace description and the instantaneous desired 
channel. The performance for this equalizer is lower than the first since the 
interference matrix includes subspace components from past interference channel 
states. The relative performance decrease is dependent on the angle of arrival spread 
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and the product of the time window length and the Doppler frequency. A third 
equalizer solution involves subspace averaging for both the desired and undesired 
channels with still lower performance. 

From the above, it would seem obvious to choose the instantaneous criteria. 
However, in time varying interference dominated channels, blind equalization 
algorithms which rely on signal property restoral are subject to false signal tracking. 
False tracking occurs when the instantaneous desired channel is in a deep fade and 
the algorithm converges to a stronger interference signal. With a time average 
interference covariance solution, the algorithm becomes less susceptible to false 
tracking as the effective interference subspace averaging window is increased. 
Thus, a trade-off exists between achieving the optimum instantaneous SINR 
performance and robust signal tracking. This motivates us to choose an SINR 
equalizer solution which has the form 

* <21) 

where (•) implies an estimated value. We note that the trade-off between 
instantaneous SINR performance and tracking robustness can easily be adjusted by 
varying the forgetting factor X. Increasing X results in a decrease in the rate of false 
tracking occurrences, and a decrease in performance as compared to the ideal 
instantaneous solution. Several techniques are available to detect and recover from 
false tracking events, so signal quality is maintained as long as the occurrence rate is 
low. The proposed equalization strategy is divided into five steps which are 

1 . Estimate the present signal sequence by passing the equalizer output 
through a decision device. [Equations (22) and (23)]. 

2. Estimate the present channel state with an optimal estimation filter. 
[Equations (25) and (27)]. 

3. Estimate the received interference plus noise signal sequence. Use this - 
error signal to update the interference covariance estimate. {Equations (30) and 
(31)]. 

4. Update the spatial equalizer weight vector given the estimated desired 
channel and interference covariance. [Equation (21)]. 
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5. Repeat the process with new desired signal estimates. 
This adaptation procedure is represented by the block diagram shown in Figure 3. 
Each of the algorithm steps are described below. 

Using the present equalizer setting, a short block of received data is 
equalized 

5 o.* = *vX = [ V^-wiJ (22) 

The equalizer output is then passed through a non-linear detection function d( ) 
which restores some known property of the transmitted portable signal. In the case 
of an FM signal set, the nonlinear function can simply force the signal modulus to a 
constant value, or a more sophisticated combination modulus forcing and phase 
filtering function has been proposed. In the case of finite alphabet digital signals, 
the nonlinear function is a symbol waveform detector. The output of the detection 
function is the estimated signal sequence 

S o,* =d < 5 o.*> (23) 

In order to accurately determine the interference subspace and the optimal 
weight vector, it is critical to have an accurate estimate of the desired channel. The 
channel exhibits highly structured time auto-correlation behavior for each antenna 
element. This structured variation should be exploited to improve our ability to 
identify and track the desired channel. The channel is estimated in two steps. First, 
a "measurement" of the reverse channel h O K is obtained by solving 

K* s **e^ n \\ h *o.r**\\ 2 r (24) 

Here, the channel is assumed time invariant over the data block. The data block 
must be a very small fraction of the inverse Doppler frequency to maintain this 
assumption. Assuming the noise plus interference is Gaussian, the solution to 
Equation (24) is the deterministic maximum likelihood estimate 

V^A',* (25) 

The channel measurements defined by Equation (25), provide a true channel plus 
noise measurement model since 
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** S o.,=V* + W *So.*=*o.» + u * (26) 

where N K =[n K — n K . N+l ] and u K is the equivalent channel measurement noise term 

The final channel estimate is a weighted sum of the present and past channel 
measurements 



n-0 



(27) 



with the FIR filter vector f given by 



(28) 



5 The solution to Equation (28) is a Wiener filter given by 

r=[£(p,p;[ + 4,p;(]-4,p;| < 29 > 

where 

is a vector of scalar fading variables representing the time varying complex 
amplitude of one channel vector entry in h K and ^ K = 

is the corresponding time varying measurement noise vector for each entry in u*. 
10 The filter solution to Equation (29) can be based on ensemble average statistics or 
based on worst case sample path behavior. The ensemble average covariance 
function for P K is: 

The noise covariance is determined by the nature of the additive receiver noise and 
the short term cross-correlation and autocorrelation behavior of the interference 
15 signals. 
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The solution for the channel estimation filter is also dependent on two 

°" 

variables, the maximum Doppler shift and the channel measurement . In the interest 
of complexity reduction, the filter design can be based on the maximum expected 
Doppler shift and the minimum expected SNR. Performance is enhanced if the 
Doppler and SNR are estimated and the filter coefficients are adapted either by 
resolving Equation (29) or by choosing the estimation filter from a fixed set of pre- 
computed possibilities. 

To estimate the undesired interference subspace, the undesired signal 
sequence is estimated 

*u,* = X ^A),* (30) 

The interference co variance estimate is then updated 

The weight vector is now given by Equation (21). The sample index k is 
incremented by N and the next iteration is initiated. 

It should be noted that many variations on this basic algorithm are possible. 
Several implementations involving sample by sample recursions have been 
evaluation. The FIR channel estimation filter can be replaced by an observer based 
filter with various state-space model possibilities. The fundamental contributions are 
the blind SINR adaption criteria, the subspace decomposition and the use of 
correlated desired channel behavior to improve estimation accuracy for the desired 
channel and the interference subspace. 

4. SIMULATIONS 

Simulations have been conducted to test the spatial equalizer approach 
discussed in the previous section. The signals are narrow-band FM with the 
following parameters: 

f c = 1GHz (RF center frequency) 
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Id = ±103Hz (v = 70mph) 

voice = 4KHz (voice waveform bandwidth) 

Afvoke = ±.12KHz (peak voice FM excursion) 

f sAT = 6KHz (SAT tone offset) 

Af sAT - Jl2KHz (peak SAT tone FM excursion) 

Here, we present one example simulation scenario in which the base station 
antenna is a six element linear array with inter-element spacing of A/2. There are 
three transmitted mobile signals, each with equal average strength. Each signal 
arrives from four different paths. Each of the four paths for each signal are spread 
over a 30° range. Each path for each signal has independent fading generated by 10 
randomly located local scatterers. The central path angles of arrival for the three 
signals are 40°, 20° and -70°. Complex Gaussian band limited noise is added at the 
output of each antenna element. The ratio of the time averaged received signal 
power to average additive noise power is 20dB per element. As the instantaneous 
channel fades, the instantaneous signal to additive noise power at each element 
varies above and below 20dB. 

The blind spatial equalization algorithm described in the previous section was 
used to recover estimates for each of the three FM signals. The detection function 
d( ) is the complex sign. The data block length NT is 200/xS with a corresponding 
adaptation rate of 5KHz. This is sufficient to support the assumption that the 
channel is constant over each block. The blind algorithm restarts when false signal 
tracking events are detected. False tracking detection can be implemented by several 
methods including monitoring SAT tone quality. 

The instantaneous SINR for the recovered signal waveforms s^ at the output 
of the spatial equalizer for each of the three signal estimates is plotted (solid) along 
with the corresponding optimum theoretical instantaneous SINR performance (dash) 
in Figures 4 through 6. The instantaneous SINR for one omni-directional antenna 
element (dash-dot) is also plotted to provide a reference for SINR improvement. 
Optimum performance is determined by the maximum eigenvalue of the 
instantaneous matrix pencil 

tV A* 
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All three simulation runs were performed on the same data set with the same initial 
weight vector setting. The only difference is the signal selected for tracking. 
Conventional CM algorithms were also applied to the simulated data set. None of 
the CM algorithms were found to converge to any of the three signals. 
5 The one second time average estimated signal SINK for signals 1 through 3 is 

18.8dB, 17.2dB and 15.7dB. The average estimated signal SINR loss with respect 
to the optimal instantaneous solution is 0.7dB, 1.9dB, and 3.8dB respectively. This 
is accounted for by the performance loss associated with the subspace forgetting 
factor A. Figure 4 illustrates the theoretical limits for instantaneous SINR 
10 performance sacrifice for signal 3 due to increasing the interference subspace 

forgetting factor A. In Figure 4, the optimal instantaneous SINR solution (solid) is 
plotted along with the maximum eigenvalue of the pencil 

(dash) with X = 0.9999 which approaches infinite averaging. The theoretical 
performance penalty, averaged over one second, between the optimal instantaneous 
15 interference covariance solution and the infinite average interference covariance 

solution is 4.6dB for the signal 3 example. These simulation results suggest that the 
proposed algorithm achieves optimal theoretical performance when the theoretical 
penalty due to interference subspace averaging is taken into account. 

CONCLUSIONS 

20 A new approach for blind adaptive spatial equalization in rapidly fading 

channels has been developed. The algorithm performs well in an extreme narrow 
band FM communication environment with high Doppler spread, high angle spread, 
and strong multiple co-channel interferers. We believe that these are the first 
reported results for a blind spatial equalizer which performs well on AMPS type 

25 signals in a high velocity portable multipath setting. The approach involves two 
fundamental results: 

1 . For certain communication signals in rapidly fading channel settings, a 
modified version of the SINR adaptive equalizer cost function is superior to 
conventional squared error cost functions. 
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2. Optimal estimation techniques can be applied to obtain accurate channel 
estimates in the presence of noise and interference. High quality channel estimates 
are crucial to obtaining an accurate estimate of the interference subspace used to 
determine the SINR equalizer solution. 

Extensions of our approach have been developed for multiple antenna delay- 
spread digital FIR communication channels. These time varying space-time 
equalizer results will be reported in a future paper. 
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WHAT IS CLAIMED IS: 

1. A base station for communicating over a radio network with a remote station, 
comprising: 

a plurality of antennas and receivers to receive inbound information signals 
from the remote station; 

a signal processor coupled to said receivers and configured to compute a 
plurality of weights for the purpose of weighting each of the receiver information 
signals and combining said weighted signals together in a manner that improves the 
signal quality at the output of said signal processor for both geographically 
stationary and rapidly moving remote stations. 

2. The base station of claim 1, wherein said signal processor includes: 
means for statistically characterizing a receive channel vector, said receive 

channel vector being representative of a receive communications channel over which 
signal energy is transmitted from said at least one remote communications station to 
said central communications station; and 

means for generating said plurality of weights based on results of statistical 
characterization of said receive channel vector. 

3. The base station of claim 2, wherein said signal processor includes: 
one or more time-varying channel estimation filters for estimating said 

receive channel vector. 

4. The base station of claim 2, wherein said means for statistically 
characterizing further include means for compiling an estimated receive channel 
covariance matrix. 

5. The base station of claim 4, wherein said means for generating said beam 
pattern weight vector includes means for averaging said estimated receive channel 
covariance matrix in order to form an estimated transmit channel covariance matrix. 
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6. The base station of claim 5, further including means for using said estimated 
transmit channel covariance matrix to form objective and constraint functions, said 
beam pattern weight vector being obtained by solving said objective and constraint 

. functions. 

7. A method of communicating over a radio network with a remote station 
using a base station having a plurality of antennas and receivers and a signal 
processor, said method comprising the steps of: 

receiving inbound information signals from the remote station via said 
plurality of antennas and receivers; 

computing a plurality of weights in said signal processor for the purpose of 
weighting each of the receiver information signals and combining said weighted 
signals together in a manner that improves the signal quality at the output of said 
signal processor for both geographically stationary and rapidly moving remote 
stations. 



15 8 - Tne method of claim 7, wherein said computing step includes the steps of: 

statistically characterizing a receive channel vector, said receive channel 

vector being representative of a receive communications channel over which signal 

energy is transmitted from said at least one remote communications station to said 

central communications station; and 
20 generating said plurality of weights based on results of statistical 

characterization of said receive channel vector. 

9. The method of claim 8, wherein said method further includes the step of: 
estimating said receive channel vector. 

10. The method of claim 8, wherein said statistically characterizing step includes 
25 the step of: 

compiling an estimated receive channel covariance matrix. 
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11. The method of claim 10, wherein said generating said beam pattern weight 
vector step includes the step of: 

averaging said estimated receive channel covariance matrix in order to form 
an estimated transmit channel covariance matrix. 

12. The method of claim 11, further comprising the step of: 

using said estimated transmit channel covariance matrix to form objective and 
constraint functions, said beam pattern weight vector being obtained by solving said 
objective and constraint functions. 
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